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We have investigated dynamic behaviors of the magnetic domain wall under 
perpendicular magnetic field pulses in ferromagnetic nanowires using micromagnetic 
simulations. It has been found that the perpendicular magnetic field pulse can 
trigger the magnetic domain wall motion, where all the field torques are kept to 
be on the plane of nanowire strip. The magnetic domain wall speed faster than 
several hundreds meters per second is predicted without the Walker breakdown for 
the perpendicular magnetic driving field stronger than 200 mT. Interestingly, the 
dynamic behavior of the moving magnetic domain wall driven by perpendicular 
magnetic field pulses is explained by charging- and discharging-like behaviors of an 
electrical RC-circuit model, where the charging and the discharging of ’’magnetic 
charges” on the nanowire planes are considered. The concept of the RC-model-like 
dynamic characteristic of the magnetic domain wall might be promising for 
spintronic functional device applications based on the magnetic domain wall motion. 
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Precise control of magnetic domain wall (DW) dynamic behaviors in ferromagnetic 
nanowires has become one of the important issues for realization of spintronic devices based 
on the DW motion[l-5|. In ferromagnetic nanowires, the DW motion generally can be 
driven by applying an external magnetic field[l, 2] or a spin current ( 3 - 5]. For instance, if an 
external held is applied along the wire axis to drive a transverse DW, a held torque from a 
cross product between the external held and the transverse DW magnetization is pointing 
out of the wire plane, generating a demagnetizing held opposite to the held torque. The 
demagnetizing held rotates the transverse DW magnetization on the wire plane, leading to 
the translational motion of the DW along the wire axisjfj. For a stronger magnetic held, 
where the held torque overcomes the demagnetizing held, an inner spin structure of a simple 
transverse DW deforms with a structural transformation into an antivortex DW structure 
9fl. This phenomenon has been well known as the Walker breakdown and the threshold held 
is called as the Walker breakdown held[l 

By suppressing the Walker breakdown, the DW speed can be made faster keeping the 
simple transverse DW structure, which implies that the DW device operation speed can 
be faster with sustaining a stable transverse DW structure. Numerous studies have been 


devoted to exploring a suppression mechanism of the Walker breakdown 
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19| . It has 


been reported that the Walker breakdown can be suppressed by using an underlayer with a 


by applying an assisting magnetic held along 
15| , by modulating the nanowire geometry with 


strong perpendicular magnetic anisotropv|ll| 
the transverse or perpendicular direction [12- 

a periodic wavy shape lfi 17], or by using a cylindrical nanowire|l8l. Il9|]. It has been known 
that the suppression of the Walker breakdown is directly related to the relative strength of 
the demagnetizing held compared to the external held torque, where both of them are along 
the perpendicular (normal) direction of hat (cylindrical) nanowires. 

In this work, we consider the DW motion triggered by the perpendicular held which 
is acting not as an assisting held additional to the external magnetic held along the wire 
axis but as an only driving field, where the held torque rotating each spin is on the hat 
nanowire plane all the time. We report our micromagnetic investigation on dynamic behavior 
of the DW motion under perpendicular held pulses with systematic variation of the hat 
nanowire width and thickness, where RC-circuit-like behavior without the Walker breakdown 
phenomenon is observed. 

We have carried out micromagnetic simulations[20] based on the Landau-Lifshitz-Gilbert 
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(LLG) equation 21]. 
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where 7 is the gyromagnetic ratio, a is the Gilbert damping coefficient, M is the magnetiza¬ 
tion, and tief / is the effective held from the total energy in the micromagnetic simulation. 

The straight ferromagnetic nanowire is set to have a length / = 2000 nm . The nanowire 
width w is varied from 20 to 100 nm and the wire thickness d is varied from 2.5 to 7.5 nm. 
The detailed geometry and dimension of the wires are schematically illustrated in the Fig. 
1(a). The external driving held pulse strength H ext is varied from 10 to 300 mT. The 
magnetic held pulse width was set to be 1 ns. Both the rise and fall times were set to be 
0.1 ns. The driving pulse held is applied along the normal direction (+z) perpendicular to 
the nanowire plane, as shown in the Fig. 1(b). In all simulations, the material parameters 
of Permalloy are used with an exchange stiffness coefficient of 13 x 1CD 12 J/m, a saturation 
magnetization Ms of 8.6 x 10 5 A/m, and zero magnetocrystalline anisotropy. The unit cell 
dimension is set to be 2.5 x 2.5 x 2.5 nm 3 and the Gilbert damping constant a = 0.01. A 
head-to-head transverse DW is initially prepared at the nanowire center before applying the 
perpendicular magnetic held pulse. 

As described in Fig. 1 , all spins (dark green arrows) are kept to be on the nanowire plane 
(xy-plane), even after the external held torque is exerted, as denoted by the dotted black 
arrows in the figure. Based on the LLG equation, it is found that the external held torque is 
t h = M x H eff , so that the magnetization direction is determined by .. y | hd x 11 „/. 

The demagnetizing held torque is 77 = M x Hdemag, so that = —|y|M x H demag (see 
Fig. 1(c)). Under a perpendicular magnetic held, inner spins of the DW are rotated by 
the external held torque along counterclockwise direction, effectively moving the transverse 
DW to the left (—x) direction along the nanowire, as shown in Fig. 1(b). Due to the 
counterclockwise rotating torque on all spins lying on the wire axis, a gradual change of 
magnetization generates magnetic source charges acting like positive ’N’ and negative ’S’ 
’’magnetic charges” [221] at the wire edges, as denoted by magenta (N) and cyan (S) colors at 
the wire edges in the Fig. 1(b). By the generation of magnetic charges, demagnetizing helds 
are generated along the —y (+y) direction on the left (right) side of the DW, as denoted in 
the Fig. 1(c). In the meantime, another torque (dotted black arrows) is generated along the 
—z direction by the demagnetizing held from the magnetic charges on the wire edges, leading 
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to the accumulation of magnetic charges on the upper and lower flat planes of the nanowire 
by the magnetization tilting toward the z axis direction, as denoted by ’N’ (magenta) and ’S’ 
(cyan) poles in Fig. 1(c). Due to the magnetic charge accumulation, another demagnetizing 
field is generated toward the — z direction as depicted in the Fig. 1(d). After the magnetic 
field pulse is applied, the role of the driving field is replaced by the demagnetizing field 
from the accumulation of magnetic charges on the wire planes, resulting in the DW motion 
toward the right (+x) direction in the nanowire, as shown in Fig. 1(d). 

To study the DW dynamic behavior driven by the perpendicular field, we have analyzed 
magnetization with respect to the time, as illustrated in the Fig. 2(a), where time-dependent 
magnetization profiles are plotted together with a temporal profile of the perpendicular mag¬ 
netic field pulse. In the figure, magnetizations are plotted for a field with 50-mT strength, 
1-ns duration time, and 0.1-ns rise/fall time. The nanowire width, w = 50 nm and the thick¬ 
ness, d = 5 nm. It is observed that the ^/-component magnetization M y does not significantly 
change under the pulse field driving since M y mainly comes from the magnetization of the 
structure-stable transverse DW, which is also evidenced by the negligible change of exchange 
energy curve observed in Fig. 2(b). The total energy is mainly dominated by the energy 
gain in the Zeeman energy when the pulse field is on. While the perpendicular pulse field is 
on, the transverse DW is simply translationally displaced along the wire with keeping the 
inner spin structure stable, as confirmed in the Fig. 2(c). Conversely, the time-dependent 
M z profile is significantly changed under the pulse field driving, temporally consistent with 
the applied field pulse, which means that the demagnetizing field is almost instantaneously 
generated toward — z direction and kept to be nearly constant for the field pulse duration, as 
also observed in the demagnetization energy curve in Fig. 2(b). It should be noted that the 
demagnetizing field generated toward — z direction is originated from the magnetic charges 
accumulation on the nanowire plane, as described in Fig. 1(d). 

M x profile exhibits a distinctive behavior compared to other magnetization components. 
Moreover, very interestingly, after the pulse field is off, the DW gradually moves back to the 
original position again with keeping the stable spin structure. As vividly demonstrated in 
the Fig. 2(a), careful analysis reveals that the time-dependent DW position profile (open 
square), determined by processing the simulated DW images (see Fig. 2(c)), is exactly 
matching with the time-dependent M x profile (red solid line). In case of Fig. 2, the maximum 
displacement of the DW under the perpendicular field pulse is 146 nmwhere the stable inner 
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spin structure is kept. . In the case of Fig. 2, the average speed of the DW motion is about 
134 m/s and could be faster than when applied held pulse becomes stronger. Under the pulse 
held with 200-rnT strength and 1-ns duration time, the average DW speed is observed to be 
about 433 m/s, where the inner structure of the transverse DW is still stably maintained. 
Reminding that the time-dependent DW position curve has the same prohle shape with the 
time-dependent M x curve, the change of the instantaneous slope from the steep one to the 
more moderate one, looks like the curve of an exponential function. If so, the DW speed is 
expected to be gradually decelerated during the pulsed time. 

To examine the relation between the time-dependent DW position prohle and the driving 
held pulse, we applied a sequential multiple held pulses perpendicular to the wire plane. 
Where 5-times held pulses with 1-ns duration time and 100-mT strength are successively 
applied. It is found that the prohle of time-dependent M x and DW position present a kind 
of continuous behavior, as shown in the Fig. 3(a). These disconnected prohles seem to be 
quite coherent and the slope change is continued even after new held pulses are applied. One 
can think that the DW motion driven by the next held pulse depends on the DW motion 
driven by the previous held pulse, suggesting that the continuous accumulation of certain 
physical quantities — ” magnetic charge”. 

Based on above confirmations and the exponential-curve-like prohle shape, we are allowed 
to assume that the phenomenon is related to a charging/discharging effect of the magnetic 
charge. The assumption is thought to be valid since there exist magnetic charges on the wire 
planes and the magnetic charges might behave in a similar way to the electrical charges. 
In order to prove our assumption, we consider the Thiele’s equation based on the LLG 
equation 23]: 

— GxX — DX + /cX = u(y x H ext ), (2) 

where, G is a gyrovector and X is a DW position. D is a damping tensor (negative value) 


related to the sample geometry, the DW width, and the Gilbert damping constant 


24]. k and 


u are a stiffness coefficient and a constant, respectively, depending on the sample geometry 
(demagnetizing factor) 251, 126]. H ex t is an external magnetic field define to be along the 
^-direction. We have found that the DW dynamic behavior is very similar to the one¬ 
dimensional (ID) version of dynamic magnetic vortex core behavior in the ferromagnetic 
nanodisk without the gyroforce term. Therefore, we try to use the ID Thiele’s equation 
—D ■ X + A;X = u (y x H ext ) to explain the observed behavior in the present study. According 
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to the ID Thiele’s equation, the time-dependent DW position function x[t] follows a relation 
as x[t\ = Ci uH ^ xt (1 — e~ kt /\ D \), where C\ is a constant. It should be noted that the resultant 
function form is very similar to the charging function form Ao[l — exp(—t/r)\ in the RC- 
circuit model j^. The role of \D\/k seems to be corresponding to an RC characteristic time 
r and the uH ^ xt corresponds to A 0 . 

In addition, we have also found that the exponential behavior of the DW position does 
not significantly depend on the rise/fall time, as shown in the Fig. 3(b), where two M x 
profiles with different rise time (0 and 0.2 ns) are compared. Both profiles are almost same 
and almost exactly matching if shifted by the difference of rise time. The change of field 
strengths does not distort the exponential behavior, either. Thus, it is concluded that the 
exponential behavior of the DW position is not originated from the field parameters such 
as rise/fall time and field strength. It is also found that the exponential behavior of the 
DW position substantially depends on the Gilbert damping constant, as demonstrated in 
the Fig. 3(c), since the RC-characteristic-time-like \D\/k term is a function depending on 
the Gilbert damping constant 24]. Every aspect of the observed exponential behavior seems 
to let us firmly believe that the RC-like DW dynamic behavior is well explained by charging 
and discharging of magnetic charges on the nanowire. 

According to the RC-circuit model, the nanowire itself can be considered as a ” magnetic 
capacitance”, where the capacitance C m is proportional to the area of the wire planes and 
inversely-proportional to the wire thickness. To check our scenario, firstly, we have analyzed 
the fully charged state corresponding to the state with the maximum DW displacement. We 
simply assume that the present ferromagnetic nanowire system is analogous to the electrical 
RC-circuit system, where an RC characteristic time constant r and the M x /Ms profile 
are fitted with an RC-circuit-like fitting function of A 0 [l — exp(—t/r)\ for charging and 
A 0 exp(—t/r ) for discharging. A 0 represents the saturated state and thus, corresponding to 
the state with the maximum DW displacement. The fitting is found to be extremely well 
matching for both the DW position and M x /Ms , as demonstrated in Fig. 3(d), where three 
fitting curves in the cases of different field pulse strengths of 20, 50, and 100 mT. Note that 
the profiles for the case of static (not pulsed) applied field become continuously exponential 
in the figure. From the fitting, A 0 values are determined to be 0.06, 0.16, and 0.33 for cases 
of 20, 50, and 100-mT perpendicular field respectively, which indicates that the magnetic 
charge at the fully charged state (Q m ) is proportional to the applied field strength (H) with a 
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coefficient of C m = ~ 3 /igA T 1 28]. Moreover, we have found that the characteristic 

times for different field strengths are r ~ 0.79 (20 mT), 0.80 (50 mT), and 0.83 ns (100 mT) 
become almost same regardless of the field strength. Considering r mainly depends on the 
wire geometry, dimension, and material parameters in the magnetic RC-circuit system, the 
negligible change of r with different field strengths is well explained. 

In Fig. 4(a)-(c), the time-dependent DW position profiles are plotted with various widths 
and thicknesses of nanowires under the magnetic field pulse of 50-rnT strength and 0.1-ns 
rise time. The widths are varied to be 20, 50, and 100 nm and the thickness is changed to 
be 2.5, 5.0, and 7.5 nm. In case of the same wire width, the DW maximum displacement 
becomes smaller for thicker wires. Quantitative analysis reveals that the wire thickness is 
inverse-proportional to the magnetic capacitance, which is explainable by the fact that the 
thicker wire has the lager separation for a plate capacitor. In a similar way, in case of 
the same wire thickness, the DW maximum displacement becomes lager for wider wires, 
where the wire width is confirmed to be proportional to the magnetic capacitance. The 
observed behavior is again explainable by the fact that the wider wire has the wider plane 
area for a plate capacitor. Therefore, the relation between the magnetic capacitance and 
the wire dimension could be simply deduced as C m oc (w xl)/d. If d, l, and w are fixed, the 
capacitance will be a constant, as C m ~ 3 /igA 1 T _1 for the case of the Fig. 3(d). 

Considering the relationship between the capacitance C m and the wire dimension (d and 
w ), we have further analyzed the maximum displacement of the DW motion under static 
perpendicular fields with variation of the field strength from 10 to 200 mT. Reminding that 
the displacement of the DW position is directly related to the generation of magnetic charges 
on the wire plane,the magnetic static field is considered to act like an energy potential as 
in the RC-circuit model. Interestingly, we have observed that the slope of the maximum 
displacement profile with respect to the static field is constant regardless of field strengths 
for nanowires with various dimensions, as plotted in the Fig. 4(d). This also matches with 
the RC characteristic time is mainly determined by R and C, not by the electrical potential. 
In magnetic analogy, the characteristic time r is also constant for each wire regardless of the 
field strengths (acting like the potential). To further confirm the RC-circuit characteristic of 
the DW motion under perpendicular pulse fields, the resistance R is calculated by the relation 
function r = RC for each wire dimension. We have found that the R value also becomes 
nearly constant for each wire, not depending on the external field strength, as shown in Fig. 



4(e), whereas, under driving field stronger than 100 rnT, the R value gradually increases, 
in particular, fora relatively narrower and thinner wire (w = 20 nm and d = 2.5 nm). This 
deviation might be explained considering that a stronger perpendicular field modifies the 
DW width, thereby changing the RC-characteristic-time-like \D\/k term which depends on 
the DW width. However, our simple RC-circuit-like model is found to be valid within a 
practical range of the field strengths not stronger than 100 mT. change, 
par 

In conclusion, RC-circuit-like DW dynamics is discovered in ferromagnetic flat nanowires 
under a perpendicular field pulse. With a systematic variation of the wire width and thick¬ 
ness as well as the field pulse strength, the RC-circuit-like dynamic DW behavior is stably 
observed without occurence of the Walker breakdown. ’’Magnetic charges” accumulated on 
the wire plane seem to well explain the observed RC-circuit-like DW dynamic behavior, 
where the displaced DW back and forth are exactly corresponding to the charging and dis¬ 
charging in the RC-circuit. The observed RC-circuit-like behavior is well explained as well 
based on a simple ID Thiele’s equation based on the LLG equation. Our study opens a new 
possibility to control the fast DW motion avoiding the Walker breakdown in ferromagnetic 
nanowires, where the DW dynamics is exactly following the well known RC-circuit-like be¬ 
havior, which is expected to be promising in designing spintronic logic devices based on the 
DW motion along ferromagnetic nanowires. 
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Figure 1. (Color online) (a) Schematics of a head-to-head transverse DW spin structure 
and a ferromagnetic flat nanowire dimension under zero field, (b) Schematics of the mag¬ 
netic charge (N and S) accumulation at the wire edges by the external field torque (dotted 
black arrows) denoted as the magenta (N) and cyan (S) colors under the perpendicular 
pulse field (H ext ). (c) Schematics of the magnetic charge (N and S) accumulation at 
the wire planes by the torque of the demagnetizing field (hollow blue arrows) along the 
transverse direction of the wire.(d) Schematics of the demagnetizing field generation along 
the -z-direction after the field pulse is off. The rotating directions (magnetization direction) 
of all spins are denoted as green arrows. 


Figure 2. (Color online) (a) Time-dependent magnetization and DW position profiles. 
H ext — 50 mT field pulse profile is on the bottom, (b) Profiles of Zeeman, demagnetiza¬ 
tion, exchange, and total magnetic energies, (c) DW positions at t — 1, 2, and 5 ns in 
the nanowire of w = 50 nm, d = 5 nm. Color code on the right represents the M y component. 


Figure 3. (Color online) (a) Time-dependent M x /Ms and DW position profiles (up) 
under multiple sequential perpendicular field pulses with H ext = 100 mT (bottom), (b) 
M x /Ms profiles with rise times of 0 and 0.2 ns. The field pulse profiles are denoted as 
dotted line, (c) M x /Ms profiles with different damping parameters of 0.01, 0.05, and 0.1. 
(d) DW positions under pulse fields (solid) and static fields (open) with 20, 50, and 100-mT 
field strengths. Fitting curves (solid line) for time-dependent M x /Ms are plotted together. 


Figure 4. (Color online) DW positions with variation of wire width and thickness. The 
thickness is varied to be 2.5, 5.0, 7.5 nm for (a) the width of 20 nm (inset is a zoomed- 
in figure corresponding to the dotted box region), (b) 50 nm, and (c) 100 nm. (d) The 
maximum DW displacement with respect to the static field strength and (e) the resistance 
R profile with variation of width and thickness, fixing the field pulse strength in each case. 
The fitting curves are denoted as solid lines. 
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